arXiv: 1509.03418v2 [cond-mat.str-el] 4 Oct 2015 


Hot-electron relaxation in metals within the Gotze-Wolfle memory function formalism 

Nabyendu Das and Navinder Singh 

Theoretical Physics Division, Physical Research Laboratory, Ahmedabad-380009, India 

(Dated: October 6, 2015) 

We consider non-equilibrium relaxation of electrons due to their coupling with phonons in a simple 
metal. In our model electrons are living at a higher temperature than that of the phonon bath, 
mimicking a non-equilibrium steady state situation. We study the relaxation of such hot electrons 
proposing a suitable generalization of the memory function formalism formulated by Gotze and 
WolfleP]. We derive analytical expressions for both dc and optical scattering rates in various 
temperature and frequency regimes. Limiting cases are in accord with the previous studies. An 
interesting feature, that the dc scattering rate at high temperatures and optical scattering rate 
at high frequencies, are independent of the temperature difference between the electrons and the 
phonons is found in this study. The present formalism forms a basis which can also be extended to 
study hot-electron relaxation in more complex situations. 

PACS numbers: 72.10.-d, 72.15.-v ,72.10.Di 


I. INTRODUCTION 

In electronic systems conductivity is one of the most 
important quantity to study. It is related to the scat¬ 
tering rate of the charge carriers due to various interac¬ 
tions. Theoretical studies of electron transport in a solid 
suffer mainly from two difficulties. One is related to the 
proper choice of a scattering rate or the inverse quasi¬ 
particle lifetime and the subtle difference in scattering 
rates found from a transport measurement, e.g. resis¬ 
tivity measurement to the one found in a spectroscopic 
measurement such as ARPESjl]. Other is to incorpo¬ 
rate the non-equilibrium nature of the charge carriers 
in certain transport measurements. In an equilibrium 
transport calculation, Drude approximation is the first 
starting point which treats electrons as non-interacting 
classical particles and an average life time or collision 
time is assigned for all. It relates the frequency depen¬ 
dent conductivity cr(a;) = croc = Here n, 

e, m and r are the density, the charge, the mass or more 
precisely the effective mass and the lifetime of the elec¬ 
trons respectively. From the expression of Drude conduc¬ 
tivity we see that the electronic lifetime or the inverse 
scattering rate and the effective mass are the most im¬ 
portant parameters to determine the transport properties 
of an electronic system. Among them, change in the r 
becomes more prominent in presence of various interac¬ 
tions and the studies of the transport properties becomes 
synonymous to the studies of the scattering rates. In the 
Drude theory, the collision time (r) is estimated as the 
ratio of average inter-particle separation(l) and the av¬ 
erage velocity(nat,g). Following equipartition theorem, it 
is assumed ~ fc^T, the thermal energy. It under¬ 
mines electron velocity and hence the scattering rates in 
a metal. Later the Fermionic nature of the electrons was 
considered and Vavg was replaced by the vp, the Fermi 
velocity by SommerfieldQ. Next comes Boltzmann equa¬ 
tion approach. In this approach, instead of looking at in¬ 


dividual particle dynamics one considers the time evolu¬ 
tion of the distribution function for the collection of parti¬ 
cles under external perturbations^. It is certainly a ma¬ 
jor improvement over the Drude-Sommerfield approach. 
However this approach is mostly limited to semi-classical 
approximations. Also the related equations governing 
the time evolution of the distribution function is mostly 
solved within relaxation time approximation. The later 
is confined to near equilibrium situations and often uses 
frequency dependent quasi-particle lifetime for better re¬ 
sults. All these traditional approaches are restricted in 
many aspects and thus fail to explain transport behavior 
in many physical systems. The time scale enters into the 
transport calculations should be the transport life time, 
which is often misunderstood as the scattering life time. 
They are certainly different. The transport properties 
are defined and hence measured in the presence of an 
external field put in a preferred direction. Thus unlike 
the case of equilibrium spectroscopic properties, scatter¬ 
ing in all the directions can not have equal contributions 
in determining transport coefficients. For example in a 
resistivity measurement in certain direction, scattering 
of charged particle in the opposite direction or the back 
scattering should have the dominating contribution. This 
difficulty can be overcome by incorporating vertex correc¬ 
tion in quantum many body theory which results in the 
famous (1 — cosdkk') contribution to the resistivity ex¬ 
pression in an electron-phonon system[2|. Here 0kk' is 
the angle between the incident and the scattered direc¬ 
tion of an electron. However since transport properties 
like conductivity, are related to the two particle correla¬ 
tions, it is never fully justified to express them in terms 
of the single particle self-energy and vertex corrections. 
The transport calculations becomes more problematic, in 
case of strong electron-electron interactions and or in the 
non-equilibrium situations. 

Here comes the importance of the memory function ap¬ 
proach to this problem. This approach was originally in- 
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troduced by MoriQ and Zwanzig separately in study¬ 
ing non-equilibrium problems. A detailed review of the 
original work can be found in Ref. Later it was suc¬ 

cessfully introduced in calculating electrical conductivity 
by Gbtze and W61fle[H. Similar approaches are also used 
recently to understand the transport properties of various 
metalic systems [lol - [l^ . They express memory function 
for electronic transport in terms of the two particle cor¬ 
relations and calculate it at the lowest order in electron 
phonon coupling. In their approach corrections beyond 
the single particle self energy calculations and the vertex 
correction, at least in the lowest order are incorporated 
naturally. 

Now let us focus on some situations where adopting 
an equilibrium description of electronic conduction where 
both the electron and phonon temperatures are same is 
not justified. In a situation where large current flows 
through a sample in presence of a large external elec¬ 
tric field or strong pulse [l3l - [l5| the current-voltage re¬ 
lation becomes non-ohomic or non linear. In such sit¬ 
uations, within an experimentally accessible time scale 
the system remains in a non-equilibrium state. This oc¬ 
curs due to the following reason. A flow of large cur¬ 
rent heats the system. Due to the lower specific heat 
of electrons than that of the phonons, the electron tem¬ 
perature becomes higher. The electrons owing to their 
shorter relaxation time scale (due to the electron-electron 
interactions) equilibrate much faster than the phonon 
subsystem. Thus they end up at a quasi-equilibrium 
state where their effective temperature is higher than 
the phonons. They are termed as “hot-electrons” here. 
They remain hot until they relax due to their coupling 
with phonons. We can also form such a quasi-equilibrium 
steady state using continuous wave laser excitations. In 
this situation electron-phonon scattering rate can not 
be calculated within an equilibrium description. One 
needs to consider “hotness” of electrons in a suitable 
way. This motivates Kaganov et. al. to propose 
the so called Two Temperature Model(TTM) which was 
later described in a modern language by Allen[i3. De¬ 
tailed discussions on this model can be found in a re¬ 
cent review by one of the authors [l^. The relaxation of 
hot-electrons due to coupling with phonons is also dis¬ 
cussed in many recent literature in a different contexts 
within different approaches [ iqI - I^ . Adaptation of such 
non-equilibrium picture got more impetus with the ad¬ 
vent of pump-probe spectroscopy. Later is recently being 
used to study important quantities like electron-phonon 
coupling in cuprates and other correlated systems ( 2 ^. 

The combined effects of these two aspects, even in a 
simple metal are not well understood. Thus we consider 
a simple metal in a non-equilibrium or quasi-equilibrium 
situation where the electron subsystem is excited to a 
higher temperature than that of the phonon bath and 
attempt a memory function generalization of the TTM 
formalism to calculate the effects of the electron-phonon 


coupling on the electrical conductivity. This is a new 
formalism which includes transport aspects in a non¬ 
equilibrium situation. We derive analytical expressions 
for both dc and optical scattering rates in various temper¬ 
ature and frequency regimes. Limiting cases are in accord 
with Bloch-Boltzmann formula as found in the previous 
studies[lj. An interesting feature, that the dc scatter¬ 
ing rate at high temperatures and optical scattering rate 
at high frequencies, are independent of the temperature 
difference between the electrons and the phonons is also 
a major outcome of our approach. We also discuss the 
possibilities of extending this approach to many other 
directions. 

This paper is organized as follows. In section [H] first 
we review the memory function Formalism used by Gotze 
and W61fle(GW) to calculated electrical conductivity in 
the subsection III Al Then in the subsection III Bl we pro¬ 
pose a new generalization of it in case of Hot-electron 
systems and derive an expression for scattering rate due 
to electron-phonon interactions. In the next section uni 
we show our results. Here we analytically calculate 
the electron-phonon scattering rate in various limiting 
cases also show it’s full behavior using numerical calcu¬ 
lations. Finally in section IIVI we summarize our results 
and present our conclusions. 


II. MEMORY FUNCTION FORMALISM 

Our aim is to study the electron relaxation in a non¬ 
equilibrium steady state situation and to calculate the 
scattering rate using the Memory function formalism. 
First we review the later in case of electron-phonon in¬ 
teraction in metals in equilibrium. Then in the next sub¬ 
section we develop a generalization of the formalism in 
case of a non-equilibrium steady state situation. 


A. Gotze and W61fle(GW) formalism[2] for 
electron-phonon interactions 


We consider a degenerate gas of electrons having 

,2 

isotropic energy dispersion Cfc = where k is the mo¬ 
mentum and m is the electronic mass. Also in the whole 
discussion we follow the system of units where h = 1 
and fcs = 1. We start with the following Hamiltonian 
describing electron-phonon interactions. 


H 

H 


Ho + Hph + H' 

+ ^ - k')cj^^Ckv6k.-k' + H. c. 

k,k' 


( 1 ) 


The hrst and the second parts of the Hamiltonian rep¬ 
resent free electrons and free phonons respectively. The 
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electron-phonon interaction is depicted by the last term, 
where the electron-phonon matrix elements is given as, 

^(q) = il/y/2miNujg)qCg. (2) 

The phonons are considered as acoustic having a disper¬ 
sion of the form ujg = cq. Here c is the sound velocity and 
q is the phonon momentum. The constant Cg = \/pf for 
metals. This completes the Hamiltonian description of 
system considered here. Before using the memory func¬ 
tion appraoch to a two temperature scenario, let us make 
a summary of the previous work by Gotze and W61fle[lj. 
In a theory of transport properties, we look for the linear 
response of a physical quantity represented by an opera¬ 
tor A due to an external perturbation coupled to another 
operator B. The response can be quantified and can be 
expressed in terms of the correlation function of these 
two operators as, 

pCO 

Xab{z) = {{A]B)) =-i I e“‘([H(t),H(0)])dt. (3) 

Jo 

Here < . > describes expectation value of the inside 

operator and << . >> describes Laplace transform 

of the expectation value of the same. Being a correlation 
function, xab{z) must maintain causality and must be 
compatible with the equation of motions obeyed by the 
operators and thus it has the following properties, 

XABiz) = {[A,B])/z for \z\^ oo. (4) 


Since the memory function is related to the correlation 
function, it has the following properties, 

M(z—T^oo) -A Ijz, AI*{z) = M{z*), 

M'{u) = 

and has a spectral decomposed form, 

M{z) = —J dioM'' {u}) / {w — z). ( 12 ) 

For vanishing impurity concentrations and electron- 
phonon coupling, the current j is nearly conserved and 
X,M -A 0 in this limit. Thus neglecting the terms 
0{xM) in Eqnirni we get, 

zx{z) =XoM{z). (13) 

Substituting A = B = ji in equation of motion (EqnIS]) 
we get, 

= (14) 

Now considering A = and B = ji we get the 

following relations, 

z((A;j,}} = ([A,j,]}-((A;A}}. (15) 


In terms of the Laplace transform of derivative, equation 
of motion is given as. 


z((A;B}} = {[A,B]} + {{[A,H];B}} 

= {[A,B])-{{A-,[B,H])). (5) 

Since xab{z) is analytic, we can write it as a spectral 
decomposition as follows. 


Xab{z) = \ j duJXABi^)/i^ - ^)- 


( 6 ) 


Our focus is on frequency dependent conductivity. It is 
related to the current-current correlation function within 
Kubo formalism!^ as, 

, ,2 


a{z) = -i—xiz)+i-^. 

Z ATTZ 


( 7 ) 


In this case correlation is considered between the com¬ 
ponents of currents in the different directions. In an 
isotropic case, the related correlation is defined as, 

{{ji'Jj}) = ji = (8) 


Now we attempt to express the frequency dependent con¬ 
ductivity in terms of a Memory function M(z). For this 
problem the later is defined by the relation 


a{z) = ^u;^^/[z + Miz)], (9) 

where ^ is the square of the plasma frequency. 

Using EqnjTl M{z) can be written as, 


M{z) 


zx{z) 

Xo - X(z)' 


Since the hrst term in the right hand side is frequency 
independent, we get in the z = 0 case, 

([Aji]) = ((^;A)).=o. (16) 

Using the above relation, we can rewrite the memory 
function as, 

XoM{z) = {({A- A)), - ((H; H)),=o) /z. (17) 

Defining (j){z) = ((H;^))^, where A = [ji,H'], the mem¬ 
ory function can be approximated in the lowest order in 
electron-phonon coupling as, 

M{z) Ki [(j){z) - (jiO)]/{zxo)- (18) 

The above expression forms the basis for calculating the 
memory function for different interactions in the GW 
original work. However our focus is on electron-phonon 
interaction and in this case M{z) can be calculated as 
follows. 

We use the expression for electrical current as j = 
SicCT''^(^)'^ko-'^k'(T, where v(k) is the velocity of the car¬ 
riers. Now from the electron-phonon interaction part of 
the Hamiltonian (EqnlT|)we get, 

A = ^[u(k)-u(k')][i:)(k-k')cj^^Ck'<7fek-k'-77-c.]. (19) 


( 10 ) 
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Using the above expression, the two particle correlator can be written in the form, 




kk' <7 pp't 

X [^[£i(k - k')][U>*(p - p')]((cLckv&k-^k'; cj,^Cpv^p-p')) + c-c. 

X ~ “ k')^[/(l - /')(1 + n) - (1 - f)fn] 


X 


1 

^k ^k' k' “I" ^ 


1 

^k ^k' ^k—k' ^ 


( 20 ) 


Here / = /(efc,T), /' = f{ek',T) are the Fermi-Dirac distribution functions at different energies and n = n{ujq,T) is 
the Bose-Einstein distribution function at temperature T. Thus we reach at the Gotze and Wolfle expression for the 
correlator(Eqn. (53a), Ref. Q), 


(j){z) 


o M ^ X ~ “ k')^[/(l - /')(1 + n) - (1 - f)f'n\ 


kk' 


Cfc ~ Cfc' 


Wk-k' 


Cfc ~ Cfc' 


Wk-k' - -2 


( 21 ) 


Using the above expression in Ean llSl followed by analytic continuation and after some algebra, the imaginary part 
of the memory function turns out to be. 


M» = |^^|D(k-k')|"(k-k')^[/(l-/')(l + n)-(l-/)/H 

X - [5{ek - Cfc' - Wk-k' + uj) - 5{ek - Cfc' - Wk-k' - w)]. (22) 

U! 


These summarizes the memory function formalism to the case of electron-phonon coupling and we will extend the 
above expression in a two temperature scenario in the next subsection. 


B. Extension of the GW formalism to a non-equilibrium steady state situation 


Now we make a generalization of the above expression to the case of TTM in a phenomenological manner. Here 
we consider the case when the temperature of the electron subsystem and the phonon subsystem are different and 
are given by Te = l//3e and T = 1//3 respectively. Thus the thermal factor before the first delta function in Eanl^ 
becomes. 


1 1 g/3eefc< g/3(efc-efc,-|-i^) ^ 1 1 

ijj -|- X -\-1 g/3(efc—Efc'-i-ii') — X uj -|- 1 -|- 1 el^i''k—^k'+‘^) — X 


(23) 


The thermal factor before the another delta function has the same structure but with w —>■ —oj. In the expression for 
the imaginary part of the memory function, the combination of the thermal factor and the delta comes as. 


/r = - Ck' - Wk-k' + w) 

UJ 

-fterms with w — ^ —uj. 


1 


1 


gPe^k _|_ X gPe^k' X e^(^k—<'k'+‘^) — 


t(< 


.P{<'k-<'k'+‘^) _ pP^if-k-f-k') 


(24) 


To simplify the Eqnl22] further, we insert an integrated delta function (1 = f dqS{q — |k —k'|)) and the resulting 
expression becomes. 


, , 2 ttN^ 

^ H = o^— 
3 TO^xo 


dTT- 


k^dk 

Stt^ 


27r- 


k'‘^dk' 

Stt^ 


dO sin0- 


1 


2miNuJk- 


Ik-k 


X J |k-k'|) X /t.(25) 


In case of metals, three of the integrals can be preformed exactly and the expression \T2\ becomes. 


M"{uj) = Mq J |X(n(/3,w,) - n(/?e,Wq - w)) (w - Wq))-f w-wj . 


(26) 
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The details of the calculations are given in Appendix!^ Here n(/3, e) is the Bose-Einstein distribution obeyed by the 
phonons living at a temperature 1//3 and having energy e. Above relation is the general result for the imaginary part 
of the memory function in a coupled electron-phonon system when electron variables are integrated out. It describes 
the scattering rate of the electrons coupled to the phonon bath when they are at two different temperatures. We will 
discuss the above expression in both the dc and frequency dependent cases in the next section. We need to perform 
the integral numerically to get the full behavior in arbitrary temperature and frequency regime. However in various 
limiting cases, analytical results can be obtained and are discussed as follows. 


III. RESULTS 


Case-I: First we consider the low frequency limit. In this case Ean]26l can be simplified as follows. 

M"{uj«u!d) = lim Mo [ | — ((n(/3, w„) — n(/?e, — w)) (w — w„))-|-a; —^ 

uj^O Jq [UJ ) 

= Mo J dqq^ -n{l3e,ujg)+ujn'{l3e,ujq)-uj'^ri"{l3e,ujq)+uj^n'"{(3e,ujg)+ 0{uj'^)) 

-\-UJ —>• —Uj} 

PQD 

= 2Mo / dqq^ (ri{/3, ujg) - n{l3e, ujg) - w,n'(/3e, ujg) - uj'^n"{l3e, ujg) - oj^utgu'"{Pe,(^q)) ■ (27) 

Jo 

The frequency derivative of the Bose-Einstein distribution, n'{ujg — uj) = and its higher derivatives appearing 

here can be written in terms of the Bose function itself and is given in the Appendix |B] In the above calculations, 
integrals can be written in terms of the dimensionless variables. It is done by scaling the momentum q by the 
temperature scales defined by Tg or T as suitable. When both the electron and phonon temperatures are lower than 
Debye temperature, the upper limit in the integral for the dimensionless variable becomes infinite. The resulting 
integrals {Dn{x)) are known as Debye integrals, as they appeared first in the specific heat calculations within the 
Debye model. They are related to the Gamma function and the Riemann Zeta function and are shown in the Appendix 
o Now we discuss the above equation in different sub cases. 

Subcase (a): We first consider a; = 0 and T,Tg << i.e when both the electron temperature and the phonon 
temperature are lower than the Debye frequency. In this case, 


M"(w = 0) = 2M, 


1 


r(5)C(5) - 


c/3)5 

= 2Mo{AT^ + BT^), 


(c/3e)' 


-r(5)C(5) 


(c/3e)' 


-r(6)C(6) 


dxx^ 


{cPef Jo - 1)^ 


(28) 


where A = ^r(5)C(5), B = -■Ar(5)C(5) -F ^r(6)C(6) -h ^ For /3g = /3, this result is identical to the GW 

result for scattering rate due to the electron phonon interaction and leads to the famously known Bloch formula for 
the resistivity in an equilibrium electron-phonon svstemj^. 

At a finite but low frequency the frequency dependent terms in Eqn. [23 contributes as follows, 

rQD 

AM''(a;-^0) = -2a;^Mo / dqq^ {rt^Pe.ujq) ^ ojqn''\Pe,^q)) 

fQD 

= -2uj‘^l3lMo / dqq^ {n{l3e,ujq) -\-3n^{Pe,ujq) -\-2n^{Pe,iOq) 

Jo 

PeOJqTl(^PeT ^q') (ojq'j \2^QUJqTl (uJq'j 3^eOJqTl 

The coefficient involves integrals with various powers of the Bose-Einstein distribution function. They can be estimated 
as follows, 


fQD 2 

Ik = dqq^Pe^^qn^i^q) = 

Jo 


X" 


{PeCf 


r/Secqa 

lo (e" - 1)'= 


dx ; 


^ ry 1 if Tg » wo¬ 


rn 


Thus, there is a contribution and the coefficient has a /3g/fe dependence which ~ Tj at the lowest order for 
T,Tg << ojo, i.e. when both the temperatures are less than the Debye temperature. The dependence shows 
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that the Landau quasi-particle nature of the conducting electrons are preserved, as expected in case of weak electron- 
phonon interaction. 

Subcase (b): In the opposite limit, namely T, Tg > wd, i.e. when both the electron temperature and the phonon 
temperature are higher than the Debye frequency, Eqn. takes the following form. 


M"(w = 0) 


2Mn 


1 iPcgoY 


{cl3)^ 4 

1 {Pecqo)^ 


D4{/3cqD) - 


1 {PeCqoY 


{cPeY 


D^dSeCqa) 


{cPef 
1 


{cPef 


4 

Pecqo 


■D4{j3f,cqD) 


dxx° 


0 


(e“ — 1)^ 


2Mo 

2Mo 


\ c/3 4 c/3e 4 5 

V 5 c/3 4 c/3e 4 


1 {PeCqD)'^ \ 
{cPef 4 ) 


(31) 


Here Hi = 2^, i?i = 4^. It is important to note that the scattering rate is independent of the electron temperature 
and is dictated by the phonon temperature only. 

If we consider the non-zero low frequency case, the frequency dependence in AM"{uj —>■ 0) is still ~ But the 
coefficient has a linear in dependence. 

Case II: For frequencies greater than Debye frequencies, i.e. u » ujd, the scattering rate given by the Eqn. [27l takes 
the following form 


M'Ylo) = Mo / dqq'^ 


_ X g-/3et^ _ X 


-I- W —>• —U! 


= -Mo 
= -Mo 
= -Mo 
= -Mo 


dh _L 

5 g-^e<^ 

' q^ 1 1 X \ 

—-' dxx^— -r 1 -I- w —?► —w, X = f 3 ujq = Pcq 


pqo ^ 


5 g /3e<^ _ X (c/3)® 


/o 


dh 1 


1 [PcqoY 


5 g _ X (c/3)® 4 

dh 1 


- 1 


Diddcqo) ) -I- w —>■ —uj 




5 g Pb^b _ X 


D4{PcqD) M- w — —uj. 


(32) 


If phonon temperature is much lower than the Debye fre¬ 
quency, i.e. for fdcqD >> Ij the above equation reduces 
to, 

(f 

-|-(a; —>• —w) 

= (f 


For phonon frequencies greater than the Debye energy, 
i.e. for Pcqn << 1, the constant part of the above ex¬ 
pression will remain same and the later part will ~ T 
instead of T®. It also notable that in both cases M"{uj) 
is independent of the electron temperature. 

Subcase c): T « oju, Te —>• oo. The frequency de¬ 
pendence will be similar to the previous case. But the 
frequency dependent part will have T® dependence. Thus 
we get various limiting expressions which are summarized 
in a table below. 
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Lo and T regime 

M"{uj) 

W = 0, Te,T <<Td 

2Mo{AT^ + BT^), A = ^r(5)C(5), B = ^ (-F(5)C(5) + r(6)C(6) + /“ dxx^^^) 

w = 0, Te, T > Td 

2Mo{Ai + B,T), Ai = !l^,B, = 

UJ » UJD, T « Td 

Mo(f+ ^r(5)C(5)) 

u! » ujd, T > Td 


UJ « UJD, T, Te < Td 

2Mo (T5r{5K{5) ^ j 

UJ « UJD, T, Te > Td 

2Mo (-b Const, x TeUj^) 

UJ « UJD, T < UJD, Te > Td 

2Mo -b Const, x TeUj^) 


However to understand the full dc and frequency de¬ 
pendent behavior, we need to perform the integral in 
Ean ]26l numerically. We choose the Debye energy ujd = 
O.OSeV = 348.14i^ to perform the integral. Now we 
present our numerical results. Possible reasons of var¬ 
ious findings will be discussed in the next sections. 



FIG. 1: Variation of the imaginary part of the memory func¬ 
tion with temperature at zero frequency at different electron 
temperatures. Here M"{0,Te,T) is scaled with M"{0,Te — 
Td,T = Td) and the temperature is scaled with the Debye 
Temperature. 

In FigH] we show the temperature dependence of the 
M''{uj —>• 0,T,Te), normalized by M"{uj —>• 0,Td,Td), 
which shows a monotonic increase with temperature. 
However there is a change in the curvature at low tem¬ 
perature as the electron temperature differs from that 
of the phonons. In Fig 121 plots of M"{uj,T,Te) nor¬ 
malized by M"{0,T,Te) at a low electron temperature 
is shown. In this case we choose phonon temperature 
to be O.OOleH ~ 12K. Here we see an increase fol¬ 
lowed by a trend of saturation with the increasing fre¬ 
quency. Also with higher electron temperature, electrons 
tend to decouple from the phonon bath and the mag¬ 
nitude of the scattering rate reduces. Figl3] shows the 
same plot at an intermediate temperature 120K. Here 
we see much less rapid variation of scattering rate with 
frequency compared to that observed in Figl2] near the 



FIG. 2: Variation of the imaginary part of the memory func¬ 
tion with frequency at two phonon temperatures and various 
different electron temperatures. Here M"{u,Te,T) is scaled 
with M”(0, Te,T) and the frequency is scaled with the Debye 
frequency. 

Debye frequency. Difference in the scattering rates also 
reduces as the phonon temperature becomes higher. Ex¬ 
cept the change in the magnitude, same trend continues 
at a temperature much higher than the Debye tempera¬ 
ture as shown in Fig 14] 

IV. SUMMARY AND CONCLUSIONS 

In Summary, we propose a generalization of the stud¬ 
ies in electronic transport by memory function formal¬ 
ism to the non-equilibrium situation. Here an electron- 
phonon coupled system is studied where the tempera¬ 
ture of the electrons and the phonons are different, as 
happens in many experimental situations. As explained 
in the previous sections, vertex corrections and the fre¬ 
quency dependence of the transport life time are also 
included in this approach naturally. Thus this approach 
includes two subtle effects, namely some non-equilibrium 
effects and the proper account of the transport life time 
or the transport scattering rates of the electronic quasi¬ 
particles. This is quite timely and much sought after in 





























FIG. 3: Variation of the imaginary part of the memory 
function with frequency at an intermediate phonon tempera¬ 
ture, 12GK and various different electron temperatures. Here 
Te,T) is scaled with M"{0, Te,T) and the frequency is 
scaled with the Debye frequency. 



FIG. 4: Variation of the imaginary part of the memory func¬ 
tion with frequency. The plot is shown at a phonon tempera¬ 
ture = 470K higher than the Debye temperature ~ 350K and 
various different electron temperatures. Here M"{u},T^,T) is 
scaled with M"{Q, Te,T) and the frequency is scaled with the 
Debye frequency. 


bye temperature which changes to the linear in T behav¬ 
ior in the high temperature. This is in accord with the 
Bloch-Boltzmaim formula and is a benchmark for our for¬ 
malism in the equilibrium limit. 2) As the electron and 
the phonon temperatures start differing, in the same fig¬ 
ure, we see a change in the curvature of the scattering 
rate below the Debye temperature. At higher electron 
temperature, due to the reduced Pauli blocking, more 
electronic states becomes available below the Fermi en¬ 
ergy and the number of scattering events enhances. As 
a result scattering rate increases. 3) As we go above the 
Debye temperature, both the phonon generation and thus 
the electron-phonon scattering gets saturated. Thus the 
effects of the higher electron temperature or the non¬ 
equilibrium effects get reduced. In this limit the dif¬ 
ference in scattering rates at different electron temper¬ 
atures becomes small. 4) In frequency dependent cases, 
we see difference in the electron and the phonon temper¬ 
ature leads to no qualitative change in the curve. But 
there is an order of magnitude change in the normalized 
values of the scattering rate at low temperature. The 
difference in magnitudes of the scattering rates at two 
different electron temperatures reduces as we go to the 
higher temperature which is also due to the saturation 
of phonon generation as explained previously. These are 
very much interesting findings from the present work. 
There are quite a few scopes of extensions on this work by 
considering, coupling with optical phonons, proper band 
structure for quasi-particles, electron-electron, electron- 
impurity interactions, external magnetic field (to study 
magneto-resistance and the non-equilibrium Hall effects), 
etc. and also in more interesting non-equilibrium steady 
state situations. These are left for future studies. This 
work can be treated as a basis for studying the effects 
of general electron-Boson coupling as well as non¬ 
equilibrium dynamics [ 2 ^ often discussed in the case of 
cuprates and other strongly correlated systems. 


context of recently flourishing non-equilibrium spectro¬ 
scopic experiments. In our results we see I) in Fig HI 
if both the electron and the phonon temperatures are 
same, the dc scattering rate has a T® rise below the De- 
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Appendix A: Derivation of Eqn. 1261 


Since we are considering metals having degenerate electrons, the Fermi energy and the magnitude of Fermi wave 
vector kp are very large. For low energy scattering, angular changes in the scattering wave vectors are dominating 
while magnitude of the scattering wave vectors can be assumed constant and ~ kp. Thus the 0-integral can be 
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considered separately and can be performed as follows, 

/•TT 

Ie= d9 sin 65{q— \/+ k''^ — 2kk' cos 6) 

Jo 

~ j d9s'm96{q — kpy/2{l — cos9j) 

Jo 

1-2 |•2kF 

= / dx6{q-kFV^)= -^dy5{q-y) = 
Jo Jo 


kp' 


(Al) 


Considering xo = Ne/ra, and inserting the above expression for Ig, we get, 

1 r°° . Cj 


2 N 

M"{UJ) = -TT 


3 Nf. {2'K)^ramik^ 


1 


OJ 


1 


nCO pOO pOO 

/ dq-^q^ / k'^dk / k''^dk' 

Jo Jo Jo 

_ g/3e(efc-e 


^Pef-k \ 1 e/3(efc— 


(5(efe - Ck' -Ulq+Lo) 


+UJ —>■ —uj} . 

Next k, k' integrals are there to be performed. We convert them as energy integrals as follows. 


(A2) 


e = 


7 2 u /•oo poo POO POO POO 

-, de = —dk = \ —dk / k'^dk / k'‘^dk' = 2m^ / ^/ede / Ve'de' ~ 2mJeF / de de'. (A3) 

2m. m V m Jg Jg Jq Jq Jq Jq 


Thus the imaginary part of the memory function can be rewritten in terms of the energy variables as, 

,oo (J2 


, .r/// N 4: N m'^ ep 
M"{uj) = -TT- 


3 Ne rrii i2n)^kp 


1 


pOO pOO pOO 

/ dq—q^ / de de' 

F Jo Jo Jo 

i 1 _ 1 / 

See _|_ X g/3ee' _|_ X g^(e-e'+w) _ X V 


,/3(e-e'+t^) _ g^e(e-e') 


(5(e - e' - ujq + uj) 


+0; —>■ —w} . 

After performing e' integral first, we get rid of the delta function and the resulting expression becomes. 


(A4) 


, e/// N 4 N rinJ ep 
M"{w) = -TT 


poo POO 

^5 / 


3 Ne rui (27r)"‘fc|, Jg 


_g/^e(e—Cijg+a;) 


. W 

+0J —>■ —w} . 


dq^q'^ 


1 


de 


1 


1 


g/See _|_ X g/Se(e—cjg+a;) _j_ X g/St^g — X 

Now we are left with e integral A, and is given as, 


' _ g/3e(<^g-<^) 


/e = / de 
Jo 


1 


1 


g/3ee -(-X/ 1 g-/3e(e-Wg+tj) _j_ X 


(A5) 


(A6) 


The lower limit of the above integral can be considered as — oo, when energy is measured with respect to the Fermi 
energy. Now we make the following substitution, 


X = dx = (dexde. 

Using the above substitution A can be rewritten and evaluated as. 


(A7) 


A = 


f-OO ^ 

dx- 


1 


xe 


/3e(—tjg+lil) 


X V ® + 1 / \ + X 


Jo Pe \X+lJ 


dx 
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■In 


X + 4 J \a;^ + l 

X + 1 


, ^ = e 


_ p/3e(l^-l^g 


X + |- 




(A8) 






















10 


Thus the imaginary part of the memory function is left with only one integral over q variable and takes the form, 


, 4 TV ep 

M"{uj) = -TT 


poq /^2 


dq^q^ 


3 TVe m* (27r)4A:| Jq w, 


(If ^ ( 

\ Lo y — 1 V 


_ ^^e{uJq-Uj) 




g/3e(^ ^q) — 2. 
_ 4 ^ AT 6f_1 


(w — ^q) I H“ W y ^ 


Using the relations Cq = l/pr, ojq = cq and defining TUq = 

pP^{u)-U)g) 


M"{uj) = Mo / 


= Mo / 





g/3e(cJ UJq) — 2^ 


(c^ — ujq'j I -|- uj —y —UJ 


1 


0^q — 2. U^) _ 2^ 


(w — ^q) I “h W —y —0 


= Mo J ((n(/3, Wq) - n(/3e, Wq - w)) (w - w,)) + w-w|. 


(A9) 


1 /g/3.K-<.) _ g/3<.A 1 ^ ^ ^ 

5 — 1 V / e/3e(i^,-l^) _ 1 ' J 


(AlO) 


Appendix B: Derivatives of the Bose function 


Various derivatives of the Bose function can be written in terms of the Bose function itself as, 


n'{ujq — uj) 


n"(u}q — oj) 

n"'(ljjq — Oj) 


dn(ujq — uj) d / 1 \ 
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0^ (l + (!>n{ujq — w) + QrP{ujq — uj)) n'{ujq — uj) 

—(dP (l + 6n{ujq — uj) + 6rP{ujq — uj)) (n(ujq — uj) + rPiujq — uj)) 

—(d^ {n{ujq — w) + 7rP(ujq — w) + 12n^{ujq — w) + 6n'^(a;g — uj)) . 


Appendix C: Debye functions 


The first Debye function is defined as. 


j-x ,n / -1 

Dl{x)= / dt^-=x-[- 
Jo e'- 1 yn 


E 


B2kX 


2k 
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for \x\ « 2 tt, n> 1, and i? 2 fc’s are Bernoulli numbers. The second Debye function is defined by. 
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for X > 0, n > 1 Sum of the above two integrals, 
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The Riemann zeta function used in the above equation is given as, 


where r(a;) is the gamma function. If x is an integer n, 

1 


OO ^,-1 

e“ - 1 


du 


C(n) = 


r(n) 


~ , 1 
du = 


— 1 


r(n) 
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fe = l 


fe = l 


(C4) 


(C5) 
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